Convergence rates for arbitrary statistical moments of random quantum circuits 
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We consider a class of random quantum circuits where at each step a gate from a universal set is 
applied to a random pair of qubits, and determine how quickly averages of arbitrary finite-degree 
polynomials in the matrix elements of the resulting unitary converge to Haar measure averages. This 
is accomplished by establishing an exact mapping between the superoperator that describes t-order 
moments on n qubits and a multilevel 5(7(4 ) Lipkin-Meshkov-Glick Hamiltonian. For arbitrary 
fixed t, we find that the spectral gap scales as 1/n in the thermodynamic limit. Our results imply 
that random quantum circuits yield an efficient implementation of e-approximate unitary t-designs. 
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Random quantum states and unitary operators are 
broadly useful across theoretical physics and applied 
mathematics. Within quantum information science 
they play a key role in tasks ranging from quantum data 
hiding 0] and quantum cryptography Q to noise estima- 
tion in open quantum systems [4H6[ . Unfortunately, gen- 
erating an ensemble of TV-dimensional unitary matrices 
which are evenly distributed according to the invariant 
Haar measure on U(N) is inefficient, in the sense that the 
number of required quantum gates grows exponentially 
with the number of qubits, n = log 2 TV So-called 
unitary t-designs provide a powerful substitute for Haar- 
distributed ensembles. Building on the notion of a state 
t-design Q , a unitary t-design is an ensemble of unitaries 
whose statistical moments up to order t equal (exactly or 
approximately) Haar- induced values @. That is, a uni- 
tary t-design faithfully simulates the Haar measure with 
respect to any test that uses at most t copies of a selected 
n-qubit unitary. Ramifications of the theory of t-designs 
§ are being uncovered in problems as different as black 
hole evaporation and fast "scrambling" of information 
efficient quantum tomography and randomized gate 
benchmarking [Io| . quantum channel capacity [ill, a nd 
the foundations of quantum statistical mechanics 12j . 

Prompted by the above advances, significant effort has 
been devoted recently to identifying efficient construc- 
tions of t-designs and characterizing their convergence 
properties 0, |a, [I3-16]. Harrow and Low established, 
in particular, the equivalence between approximate 2- 
designs and random quantum circuits as introduced in 
[3], and conjectured that a random circuit consisting of 
/c=poly(n, t) gates from a two-qubit universal gate set 
yields an approximate t-design [16j. While supporting 
numerical evidence was gathered in (l7| for low-order 
moments, and efficient constructions of t-designs were re- 
ported in [l6| for any t = 0(n/ log ri), the extent to which 
random quantum circuits could be used to implement an 
approximate t-design for arbitrary, fixed t remained open. 

In this Letter, we address this question by determining 
the rate at which, for sufficiently large circuit depth, sta- 
tistical moments of arbitrary order converge to their lim- 
iting Haar values. Our strategy involves two steps: first, 



for given t, we show that the asymptotic convergence rate 
is determined by the spectral gap of a certain superoper- 
ator, which encapsulates moments up to order t; next, we 
compute this gap by mapping the t-momcnt superoper- 
ator to a multilevel version of the Lipkin-Meshkov-Glick 
(LMG) model, which is known to be exactly solvable, 
and whose low-energy spectrum is well understood in the 
thermodynamic limit n — > oo [3] ■ Our approach ties to- 
gether t-design theory with established mean-field tech- 
niques from many-body physics, extending earlier results 
by Znidaric 1A\ for t = 2. Furthermore, asymptotic con- 
vergence rates allow us to upper bound the convergence 
time (minimum circuit length, k c ) needed for a desired 
accuracy e relative to the Haar measure to be reached. 
For any fixed t, we find that the scaling k c ~ nlog(l/e) 
holds for sufficiently large n and small e. 

Moment superoperator.— Let a random quantum cir- 
cuit of length k be a sequence f/fc . . . U\ of k unitary op- 
erators on an n-qubit Hilbert space = ®^ l K qj , where 
each Ui is selected from an ensemble {fj,(U),U}, for a 
probability distribution /i with support on a universal 
gate set. To analyze arbitrary t-order moments, we in- 
troduce a Hilbert space "Kut = K® 2t , which consists of 
2t copies of "K and we refer to as the moment space, with 
dim(J£M t ) = D = N 2t , and a local moment space !K; t , 
which results from grouping factors corresponding to the 
same qubit in "K Mt . That is, "K Mt = ®"?fg 2t = 'KfJ 1 , 
with dim(!J-t; t ) =(1 = 4*. Moments of order t may be de- 
scribed in terms of the following linear operator on JCm 4 : 



MM 



d^{U)U fl 



Physically, M 4 [^t] may be viewed as the superoperator 
induced by the action of t copies of the random circuit 
on nt-qubit density operators defined on IK®* (see Fig. 
1). In line with standard practice in open-system the- 
ory [3] , we shall introduce "operator kets" in 'Kut i de- 
noted | A}} = A, and, correspondingly, ((A\ = A^ . Thus, 
a D 2 -dimensional operator ket transforms according to 
U®t,t\A)) = UAU\ under U £ U(N). Once a basis for 
3iu t is chosen, the matrix representation of spec- 
ifies a complete set of t-order moments. 
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FIG. 1: The moment space CKm 4 may be visualized as an ar- 
ray of 2nt qubits, in such a way that t copies support a ket 
in the state space on nt qubits, and the remaining t copies 
the corresponding bra. In this way, a unitary U on n qubits 
induces a transformation U mj on density operators on nt 
qubits. Dashed rectangles indicate the 2t qubits correspond- 
ing to a local moment space M; t , whereas the ovals correspond 
to a unitary U acting non-trivially on the first two qubits. 



The probability distribution that describes a random 
circuit of length k, is given by the k-th fold 

convolution of (i with itself [l6|. That is, ^k(U) = 
SUi=i dv(Ui)5{U - H!? =1 Ui). ft then follows that 

mm = J nUdm^Uu?^ 

= II / Mu^u®^ = (MMf = m*H 

Note that, under the assumption that the ensemble 
{fi(U),U} is invariant under Hcrmitian conjugation, 
fi(U) = n(U'), Mt[fj] is an Hcrmitian operator on "Km^ 
If n(U) has support on a universal set of gates, then the 
measure over the random circuit converges to the Haar 
measure on U(N) in the limit of infinite circuit length 
H, lim^oo Hk(U) = ^h(U). We begin by characterizing 
how these convergence properties translate in terms of t- 
order moments. Let M t [/i#] = J d/ijj(J7)J7®*' 4 , and let 

V t = span{|0)) 6 X Mt | = \tf>)), VU E U(N)} 

be the subspace of fixed points of U e U(N), with 

Tv t denoting the corresponding projector. We claim that 



lim (M t M) =9v t =M t \iHi}, Vt. 



(1) 



While this is implied by the results in [T(|, a self- 
contained proof follows. Let \<f>)) be an eigenoperator 
of Af t [/x] with eigenvalue A, and \cj>u)) = U®*'*^)). Since 



\{{MM))\ 



it follows that |A| < 1, with equality holding if and only 
if t/® 4 '*^)) = for all U with fj,(U) ^ 0. Any such op- 
erator ket \<j>)} is also invariant under any unitary of the 
form J/®*' 4 , where U is generated by a random circuit of 
arbitrary length, that is, U = H^ =1 Ui for any k, as long 
as n(Ui) 7^ 0. Thus, if fi(U) has support on a universal 
gate set, the eigenspace of eigenvalue I is precisely V t . 



Since all other eigenvalues of Af t [/i] have magnitude less 
than I, Af t [/i] converges to 7y t . To establish the sec- 
ond equality in Eq. (TT]), we invoke the invariance of the 
Haar measure under U(N), Hh{U) = pL H (U'U). For \<f>)) 
an eigenoperator of Aft [/iff] with eigenvalue A, it follows 
that, M t \ji H ]\4>)) = f dnHWU^m = A|0)). Thus, 

u^xm = Jd^iuwur^)) 

= fd i i H {U^U)U^ t \<j>)) = (dii H {U)U^\<j>)) = \\<t>)). 



If A 7^ 0, it follows that \<j>)) <G Vt, otherwise A = 0, which 
establishes the desired result. 

Our next goal is to obtain the rate at which M^[fj] 
approaches Aft Since Aft [/iff] projects onto the 

eigenspace of Af t [/i] of eigenvalue 1, the distance 
||Af t fc [/i] — Aft|/tH]|| with respect to any (unitarily in- 
variant) norm depends only on the remaining eigen- 
values {Xi} of Af t [jU] and the corresponding eigenpro- 
jectors {n^}. Specifically, if k is sufficiently large, 



\M^)-M t [^ H }\\ 



A„ fc n 



I Ai | fe ||Hx|| , where 



Ai = I — A t is the subdominant eigenvalue of Afj[/i]. 
Thus, the asymptotic convergence rate is entirely deter- 
mined by the spectral gap At of Aft [fi\ . 

Mean-field solution.— The starting point for mapping 
Aft[/x] to an exactly solvable, infinitely-coordinated model 
is to ensure that the following conditions are obeyed: (i) 
The applied quantum gates consist only of single- and 
two- qubit gates selected according to a distribution fl(U) 
on Cf(4), with fl(U) = (ii) The target pair of 

qubits is picked uniformly at random. We shall gener- 
ally refer to the class of circuits obeying (i)- (ii) as per- 
mutationally invariant random quantum circuits. Since, 
in each application of a random gate U to a fixed pair 
of qubits, the operator acts non-trivially only on 

the associated bi- local moment space 0ii t <£> r Ki t , and this 
action is identical for every qubit pair, the moment super- 
operator for any such circuit may be written as follows: 



MM 



n(n — I) 



i<j=l 



(2) 



where for any pair i,j the restriction m t [p] of Wt J [/2] 
to Wi t ® "K h acts as m t [p\ = Jd/t(C/)f/®*'*. Recall- 
ing that dim(J{j t ) = d, M t [/j] thus defines a qudit 
Hamiltonian, which is invariant under the symmet- 
ric group S„ of permutations of the n local moment 



spaces. Explicitly, if {b l a . 



*))((/? | } denotes an 



outer-product basis for operators acting on any ?£; t , 
we may expand m\ 3 = E a /3 7 a=i(( a 7l^l/3£))&a/3^,5 = 
E^ 7 fci CaPysbltpfyg, and rewrite M t [n] as a quadratic 
function of the collective operators B a p = Y^h=i ^a8) that 

is, Aff[/i] = » ( „ 1 _ 1) Yf a p 1 S=l C aPls(. B apB 1 s - 5faB a s). 

Since the operators B a p obey SU (d) commutation rules, 
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[B a p,By S ] = B aS 6^ - BfoSap, "K Mt carries the (re- 
ducible) collective n-fold tensor product representation 
of SU(d), and M t \p\ provides a d-level extension of the 
standard, spin- 1/2 LMG model (20|. 

Thanks to the invariance under S„, each of the eigen- 
operators of M t [fi] belongs to an irreducible representa- 
tion (irrep) of SU(d). Our first step is to show that the 
cigenspace V t of M t [fi] corresponding to the ground-state 
(extremal) eigenvalue of 1 lies in the totally symmetric 
irrep, of dimension d s = ( 4 * + I "" 1 ) & Recall that V t 
consists of operators in !Kjvf t that commute with all t- 
fold tensor power unitaries ?7®'. By Schur-Weyl dual- 
ity [TIED], every such operator is a linear combinations 
of elements of § t , under the natural representation in 
5£®*. Note that the operators spanning V t are permuta- 
tions of the t copies of J£ rather than permutations of the 
n local moment spaces J£j f . One may write any such per- 
mutation as |cr(™))) = Ei 1 ...-i t =iHl---^)(V(i)---V(t)|, 
where a G S t . Furthermore, each such permutation may 
be viewed as a product ket relative to the factorization 
% Mt = %f t n . Explicitly, = fla)))®", where 

W)) = Y., ..... Ni ■■■it) (V(i) ■•• I e 5Cj t . The fact 
that V t is exactly spanned by product states is significant 
from the perspective of mean-field theory. For arbitrary 
SU(d) quadratic Hamiltonians, it has been rigorously es- 
tablished that the exact ground-state energy is given in 
the thermodynamic limit by a mean-field Ansatz equiva- 
lent to assuming that the ground state is an SU (d) coher- 
ent state [20( | . Since, for the completely symmetric irrep, 
the manifold of coherent states consists precisely of all 
product states [22|, the mean-field extremal eigenspace 
of M t [/i] is, in fact, exact for any n. 

The next step is to determine the lowest excitation 
energy in the large- n limit, which is accomplished by 
expanding M t [fj] around an arbitrary extremal mean- 
field state for each irrep [HI- While, to our knowl- 
edge, a rigorous justification of such a mean-field Ansatz 
is lacking, its validity for LMG Hamiltonians is sup- 
ported by an extensive body of theoretical and numer- 
ical investigations [l8j]. For the totally symmetric irrep, 
the required diagonalization procedure is most straight- 
forwardly carried out by realizing the U(d) algebra in 
terms of d canonical Schwinger boson operators {a Q , aj^} 
[23j . That is, we let B a p — a^ap and rewrite the LMG 
Hamiltonian as M t [/x] = n ^_^ Ea/3 7 5=i c^sa^a^apas. 
Since the totally symmetric irrep of 'Kut contains ex- 
actly n Schwinger bosons, it is possible to eliminate 
one boson mode by regarding it as "frozen" in the vac- 
uum for a generalized Holstein-Primakoff transformation 
p3| . Specifically, let the local basis be chosen so that 
the frozen mode corresponds to \a)), and let 9{n) = 
(« - ala a ) 1/2 , with a\ -> 6{n), a a -> 6{n). Two 

simplifications may now be invoked: first, the fact that 
\cr' n ')) is an exact ground state causes any coefficient 



gates) to vanish; second, only terms up to the lead- 
ing order in 1/n need to be kept in 6(n). This finally 

yields: M t [fi] = 1 - ± Y? a p=i E a pa\flp + 0(l/n 2 ), where 
E a p = 2(S a/ 3- ((aa\m t \a/3)) - ((aa\m t \/3o-))). To leading 
order, the desired gap is then determined by the smallest 
eigenvalue, eti, of E a p. That the latter is nonzero may be 
shown by exploiting basic properties of the superoperator 
mt [Hi- This establishes our first main result: For any 
permutationally invariant random quantum circuit, and 
for any fixed t > 0, the spectral gap may be expanded as 



A t = ^ a p n 



n \n z 



(3) 



of the form c Q 



(and their complex conju- 



for coefficients {a p } that may in general depend on t. 

A stronger result may be obtained for a sub-class of 
random quantum circuits which arc, in addition, locally 
invariant, that is, (m(U) is invariant under the subgroup 
U(2) x U(2) C U(4) of local unitary transformations on 
the two target qubits. In this case, it is possible to choose 
a basis for each local moment space J£j t , which includes 
a maximal set of £-qubit operators in the commu- 

tant of Z7®', with U G U(2). Accordingly, every matrix 
clement ((af3\m t \j6)) = 0, unless each local basis cle- 
ment is itself an invariant, and the largc-n behavior of 
the gap is determined by matrix elements of the form 
((au\mt\crui)) and ((o~u}\mt\uJO~)) , with a G St (without 
loss of generality, we may choose \u)) = |/))) and \u)) 
an arbitrary C/(2)-invariant with ((a\uj)) = 0. Since, for 
t > 1, the maximum value of any such matrix element is 
independent of t (see [24[ for full detail), it follows that 
the leading order term a\ does not depend on t for locally 
invariant random quantum circuits. 

Example.— Consider the simplest case where t = 2 
and p.(U) = /j,h(U) on U(A). The invariant eigenspace 
V 2 of M 2 is spanned by the identity |/( n ))) = (|/)))® n 

and the permutation |5( n ))) = (\S))) that swaps the 
t = 2 copies of Jf = "K® n . Since jj,(U) is the Haar mea- 
sure, m 2 coincides with the projector onto the subspace 
V 2 for n = 2 qubits. An orthogonal basis for V 2 may 
be formed by taking even/odd linear combinations under 
swap, \A±)) = \I^)) ± l^ 2 ))). To find the excitation en- 
ergies, we choose one of the extremal local kets, |7)), and 
minimize i^inin = 2mm(l— {{Ia\rri2\Io()) — {{Ia\rri2\cxl))), 
over all local operators \a)) G "Ki t orthogonal to |7)). 
This yields |«» = \S))-((S\I))\I)) = <j\al+a\al+a\al 
and A t = 6/5n + 0(l/n 2 ). To determine how quickly the 
large-n scaling sets in, the fully symmetric sector of A7 t 
under §„ was numerically diagonalized. Since Hh(U) is 
invariant under U(2) x U(2) transformations, !K; 2 may 
be restricted to the subspace of SU (2) invariants. From 
angular momentum theory [2l| . the number of such in- 
variants is j m 2 j — (t+rjw. = ^t> where mj is the mul- 
tiplicity of the SU (2)-irrep with total angular momentum 
J. This yields d l § c = ( Ct+ „™ _1 ) < d s , which makes nu- 
merical comparisons tractable for small t. Exact results 
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FIG. 2: (Color online) Inverse spectral gap A^ 1 of M t [fijj] 
with t = 2, 3 for a random circuit consisting of two-qubit 
gates selected according to the Haar measure on £7(4). The 
line with slope 5/6 corresponds to the asymptotic result. 



for t=2 and 3 (see Fig. [2]) indicate that the scaling pre- 
diction for At becomes very accurate for n > 14. 

Convergence time.— In order to establish the use- 
fulness of a random circuit as an e-approximate uni- 
tary t-design 0, [3], we need to upper-bound the cir- 
cuit length required to achieve a specified accuracy, 
e. Let the convergence time with respect to a given 
norm be defined by the minimum length k c for which 
\\M^[fi] - M t [fjt H ]\\ < e. That M t [p] be operationally 
indistinguishable from M t [pu] requires that the supre- 
mum of ||(M t fe [/i] - M t [p H ]){p)\\i be sufficiently small 
over all nt-qubit density operators p. We may bound 
the 1-norm starting from the 2- norm (lit ]. For any den- 
sity matrix p, ||(M t fe [/x] - M t [p H })(p)h < This fol- 
lows from normalization of p and the fact that Mt{pH\ 
projects onto the eigenspace of eigenvalue 1 of M t [p]. 
In conjunction with the Cauchy-Schwartz inequality, this 
implies || (M t fe [/z] - M t [p H ])(p)\\ l < 2 nt \\. Requiring that 
2 nt \ k ^ < e finally yields k c < A^ 1 (log(l/e) + nt log(2)). 
Since, using Eq. ([3]), A^ 1 = Y^Li a ' P n2 ~ P ~ to 
leading order, k c = a^ 1 nlog(l/e) for sufficiently small e. 
It is worth stressing that we have not addressed how k c 
scales at fixed n while letting t — > oo. To answer this 
question it is necessary to determine how a[ depends on 
t, as well as the minimum value of n required for the lin- 
ear term a! x n to dominate. While this is important for a 
full characterization of t-designs, our results are directly 
relevant to physical applications, where t is fixed. 

In summary, we have established that a large class of 
random quantum circuits are efficient e-approximate uni- 
tary t-designs for arbitrary finite t. The fact that the 
extremal eigenoperators are separable suggests that sim- 
ilar results might be established for more general random 
circuits for which the Hermiticity and the §„-invariance 
assumptions of the moment superoperator need not hold 
[25I ]. As mentioned, a remaining open question is to de- 
termine how the circuit length scales as the limits of large 
n and large t are taken together. This may resolve the ap- 
parent paradox that while Haar random unitaries are in- 



efficient, arbitrary t-designs are not, possibly with equal 
asymptotic rates. We believe that our findings, along 
with the techniques introduced to analyze moments of 
unitary ensembles, will enable further understanding and 
applications of t-designs across quantum physics. 
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Appendix: Supplementary Material 

Additional remarks on the determination of the 
spectral gap 

The diagonalization procedure illustrated in the main 
text determines both the extremal and neighboring eigen- 
values of M t belonging to the totally symmetric irrep of 
SU(d). Although expressed in terms of bosonic opera- 
tors, the procedure is equivalent to a variational Ansatz 
whereby the trial wavefunction of the extremal state is 
of the form 

|0))® n = |0-*». 

and the first excited state, corresponding to a single 
bosonic excitation, is of the form 

|n* = n-l,na=l» = A=(\<l> ■■• <W) + ■■ - + W ■•■ <£»)• 

In principle, it is possible that the subdominant eigen- 
value may lie instead in the SU(d) irrep which carries 
exactly one anti-symmetric pair of indexes (see e.g. Ref. 
[21]). This can be accommodated by using a different 
variational Ansatz for the excited state to be minimized. 
If, as in the main text, we choose a local basis in J£^ t 
that includes the fixed extremal permutation | cr) ) , with 
the remaining local basis operators \a)) 7^ | cr) ) treated 
as excitation modes, the relevant single-excitation band 
is spanned by kets of the following form: 

\n a = n- l,n a =l)) = -j=(\a...o))\aa)) - \aa))), 

with nonzero matrix elements: 

«n Q = l\M t \n p = 1)) = 8 a[j - E a0 /n + 0{l/n 2 ), 
E a p = 2(1- ((aa\m t \<j(3))). 

(Note that the "exchange term" ((aa\mt\(3a)) is no 
longer present). Upon diagonalization of E a p, the sub- 
dominant eigenvalue a\ is determined, with an identical 
1/n scaling as found for the symmetric irrep. The possi- 
bility that the subdominant eigenvalue lies in this irrep 
may be removed a priori by imposing a natural additional 
restriction on the random quantum circuit, namely by re- 
quiring that the two-qubit gate distribution be invariant 
under the transformation S that swaps the two-qubits, 
that is, fi(SU) = fi{US) = fi(U). 

From a physical standpoint, it is also interesting to 
note that the existence of degenerate mean-field ground 
states indicates that M t [/i] describes a deformed (or 
broken-symmetry) phase. For generic SU(d) models in 
the broken phase, the degeneracy of the ground state 
manifold is known to be lifted at finite n by tunneling 
between the mean-field ground states, which induces cor- 
rections of order exp(— an), and hence a gap that closes 



exponentially in the thermodynamic limit, sec e.g., C. 
M. Newman and L. S. Schulman, J. Math. Phys. 18, 23 
(1977); S. Dusel and J. Vidal, Phys. Rev. B 71, 224420 
(2005) [Ref. [18] in the main text provides additional 
relevant literature]. Remarkably, no such correction can 
occur in our case because each extremal mean-field state 
|cr ( -"- ) )) is an exact eigenoperator of Mt[/x] for arbitrary 
finite n, as stressed in the text. 

Proof that the leading-order coefficient is 
non- vanishing 

We show here that, for the class of arbitrary permuta- 
tionally invariant random quantum circuits, the leading- 
order coefficient a\ of the 1/n term in the spectral gap 
expansion is nonzero. 

Recall that a\ is given by the minimum eigenvalue of 

E ap = 2(1 - ((aa\m t \a/3)) - {{aa\m t \pa))), 

where \a)) and |/3)) are operators orthogonal to a permu- 
tation \o~)) £ St. Thus, showing that a± > is equivalent 
to showing that ((o-a\m t \o~a)) + {{aa\m t \aa)) < 1, for 
any operator \a)) such that ((a\o-)) = 0. Let 

\4,)) = -j=(\va)) + \a*))). 

Upon taking the expectation value with respect to mt, 
we have 

((ip\m t \ip)) = \ f{{aa\m t \aa)) + {{aa\m t \aa}) 
+ ((ao-\m t \o-a)) + ((aa\m t \aa))^j . 

Since M t (and hence mt) is invariant under interchange 
of any two qubits, ((acr\mt\cra)) = ((o~a\mt\ao~}) and 
((aa\mt\ao-)) = ((aa\m t \aa)) , yielding 

((■0|mt|-0)) = ((aa\mt\o-a)) + ((aa\mt\aa)) . 

Following from the properties of the invariant subspace 
V t established in the main text, {{ip\mt\i}})) = 1 if and 
only if the operator is invariant under arbitrary uni- 
tary transformations of the form U® 1 with U £ U(4), 
given that the corresponding two-qubit gate distribution 
ft,(U) is universal on U (4). Thus, we must show that there 
exists a unitary transformation U £ U(4) such that 

U® t {a <Eia + a(E) er)£/ t ® t ^atgia + afgia. 

(Recall that, in our notation, we identify \a)) = a, and 
so one). In fact, we may take the permutation \o~)) to be 
the identity |/)) without loss of generality. This follows 
upon noting that 

«V>K|V)) = / dfiiUM^u^tpu^]. 
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Since aU m = U® l a for all U & 17(4), tfaU^WipW®* = 
^,t[/®*-0[/t®*. Thus, cr <8 a + a ® cr has the same expec- 
tation value as I <8 5 + 5 <8 7, where a = a^a. 

Any operator |5)) orthogonal to the identity may be 
expanded a = J2v c v< J l l ■■■ a t t i where v = (vi, . . . , v t ), 
and i>i G {0,1,2,3}, with oq = I, u\ = a x , 02 = cr y , 
(T3 = <7 Zl and the sum ranges over all possible strings v 
except the one where v% = Wi. Thus, 5<87 = J2v c^ "^ ® 
I 1 . . . er* t £8 7*. Now, under (7 = exp(|i£7 3 ® 03), the 
following transformations hold for Pauli operators: 

ffl® / H> -cr 2 (8) (T3, 

(72 ® 7 I ^ CTl ® £73 , 

(T3 <8 7 1— > CT3 ® 7, 

7 (8) £71 h-» -£73® £72, 

7 (g) £7 2 l-> C3 ® CTl , 

7 ® £7 3 I ^ 7 <g) £7 3 . 

Thus, for any |5)) with support on a Pauli string such 
that Vi = 1 or 2 for some i, there are terms in the expan- 
sion of U® t {I®a+a®I)W® t which contain factors of the 
form £7 2 (8> £73 , £73 (8 £7 2 , £7i (8 £73 , and £73 (8 £7i . Since there arc 
no such term in the expansion of 7 (8 a + 5 (8 7, it follows 
that 7(85+5(87 is not invariant under (exp(^ic73<g>£73))®*. 
Similarly, if |5)) has support on any Pauli string such that 
Vi = 3 or 2 for some i, then 7(85 + 5(87 is not invariant 
under (exp( \io\ ® 0i))®*. Since |5)) must belong to one 
of these two cases, the proof is complete. □ 

Proof that the leading-order coefficient is 
t-independent 

We next show that, under the additional assumption 
that p,(U) is invariant under the subgroup (7(2) x (7(2) C 
17(4) of local unitary transformations on the two target 
qubits, the leading order coefficient a± of the 1/n expan- 
sion does not depend on t. 

Recall that a\ is determined by the maximum value of 
((aw I mt\ cw)) + {{<juj\mt\uja)) , where a <E St and \uj)) is 
a (7(2) invariant which is orthogonal to \cr)). As shown 
previously, we may take \a)) to be the identity without 
loss of generality. Now, any operator orthogonal to the 
identity may be expanded as u; = c v o\ . . . £7* where, 
as before, v = [y\ y . . . , i>t), and vi € {0,1,2,3}, with 
(To = 7, £7i = <r x , (72 = <Ty, cr 3 = <t z , and the sum ranges 
over all possible strings v except the one where v\ = Vi. 

The first step is to show that the expectation value 
((7w|(7® t '*|7w)) for any U G (7(4) may be written as a 
symmetric polynomial of the form 

((7w|(7® t ' t |7^» = wpx Pl y P2 z p \ (4) 

2 <Pl+P2+P3<* 

where x, y, z are real numbers in [—1, 1] which depend 
only on (7, and wp are positive coefficients which de- 
pend only on with p = {pi,P2,P3) and pi being 
non-negative integers subject to 2 < p\ + P2 + pz < t. 



To establish Eq. §4§ above, we exploit the fact that an 
arbitrary element of U (4) may be written in the canon- 
ical decomposition form U = U\ ® U2 U(q, r, s) U[ (8 (7^, 
where U\, U[ , (7 2 , U'2 act locally on cither of two qubits 
and (7(g, r, s) = exp{i(£/CTi®CTi+r£72®£72+S£73!8>£73)} [see 
e.g. B. Kraus and J. I. Cirac, Phys. Rev. A 63, 062309 
(2001)]. Since \u)) is a (7(2) invariant, it then suffices to 
consider the action of U (q, r, s). Direct calculation shows 
that under U(q,r,s), the following transformations are 
obeyed by Pauli operators: 

7 (8 £7i H> riSi 7 (8 (7i + r 2 si £72 <8 £73 

- ris 2 (7 3 (8 £7 2 + r 2 s 2 <T\ <8 7, (5) 

7 (8 £7 2 ' y SiQi 7 (8 £72 + «2£7l (73 (8 CTl 

- Sl02 0*1 ® £7 3 + S 2 £72 0"2 <8 7, (6) 
7 (8 £7 3 I ^ Qin 7 ® £7 3 + ^2^1 CTl (8 (72 

- £217-2 £7 2 ® £7i + £72^2 £73 (8 7, (7) 

where gi = cos(2g), (72 = sin(2q), and similar expressions 
hold for ?'i,r2, and si,S2, respectively 

The idea is now to evaluate U m,t \Iu>)) term by term 
in the expansion for \ui)), that is, we evaluate 

l7®*(7 (8 u) U m = J2 c y U® 1 (I 1 (8 < ... 7* (8 < ) U m 

t 

= ^c,(g)C7(r®4 i )r7t, 

where now the transformation rules in Eq. ([5 ]) -([7 | may be 
applied to each of the t factors independently Comput- 
ing the matrix element ((7w|(7® t '*|7w)), there is a contri- 
bution of the form | c v \ 2 x Pl y P2 z P3 , arising from each of the 
terms c l/ 7 1 (8 ct* ... 7* (8 £7* , where x = r\S\, y = Siq%, 
z = qiT\, and pi, P2, Pz are the number of instances 
where Ui = 1,2,3, respectively. Summing over all terms 
in the expansion for \u>)) finally results in a polynomial 
of the form stipulated in Eq. (j4]), with wp = Yl v \°v\ 2 
determined by the sum over all strings v that share the 
same p- vector. That the polynomial is symmetric under 
the interchange of x, y, and z follows from the invariance 
of \ui)) under (7(2). Note that by construction, x, y, and 
z are bounded between [—1, 1]. 

Let the degree of a Pauli string be the number of in- 
stances where Ui 7^ 0. Since, under U m with U £ (7(2), 
a Pauli string can only be mapped to a Pauli string of 
equal degree, it follows that any (7(2) invariant whose 
expansion contains terms of differing degree can be writ- 
ten as a linear combination of U (2) invariants each con- 
taining only terms of equal degree. For t = 1, the only 
(7(2) invariant is the identity. For t = 2, there is exactly 
one (7(2) invariant orthogonal to the identity, namely, 
u = -^(a\af + a\a\ + o\a\). Consequently, no U (2) in- 
variant contains a term of degree 1, and the only degree-2 
terms a (7(2) invariant may contain are linear combina- 
tions of the form cr\(j\ + £720"2 + (J \ <J \- Thus, for a mono- 
mial occurring in ^ ~wpx Pl y P2 z P3 , 2 < p\ + P2 + P2 <t, 
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and if pi + j?2 + P3 = 2, then exactly one of pi, P2, or 
P3 = 2. 

The next step to establish the claimed result is to 

show that J2p w P xPl y P2zP3 ^ \( x2 + V 2 + z 2 ), where 
the right hand side is the polynomial corresponding to 
((/cj 2 |^® M |^2)), where |w 2 » is any degree-2 U(2) in- 
variant. To show this we first show that the average over 
each set of monomials x pi y P2 z P3 defined by a set of inte- 
gers p > p' > p" distributed in every distinct way to pi , 
P2 , and P3 , is less than or equal to ^(x 2 + y 2 + z 2 ) . There 
are two cases to consider: 

(i) If P > 2. then from |x|, \y\, \z\ < 1, it follows that 
x p \(y p 'z p " + y p "z p ') < x 2 , yP\{zP' X P" + Z p" ' X p ') < y 2 , 
and z p \(x p y p + x p y p ) < z 2 . Thus, the average of 
the left hand side of each inequality, which is the average 
over the desired set of monomials, must be less than or 
equal to the average of the right hand sides, which is 
\\x 2 +y 2 + z 2 ). 

(ii) If p = p' = p" = 1, then using x = fisx, y = Sx5i 
and z = qir±, xyz < \{x 2 + y 2 + z 2 ) can be writ- 
ten g 2 r 2 s 2 < \(r\s\ + s 2 q 2 + q 2 r 2 ), which follows from 
q\r\s\ < r 2 sl,s 2 qf,q 2 rl. 



Since J^p- wpx Pl y P2 z P3 is a weighted average of the above 
monomial averages, each of which is less than or equal 
to \(x 2 + y 2 + z 2 ), it follows that J^p w p xPi V P2 zPa < 
\{x 2 + y 2 + z 2 ). 

Finally, the steps described above can be applied 
to the exchange term ((TwlC/®*'*^/)), resulting in 
((Ja>|tf®*''|wl)) < \{x 2 + y 2 + z 2 ), where x = b 2 c 2 , 
y = a2C2, z = a-ibi for any U € U{4). Since these in- 
equalities hold for every U e U(4), it follows that 

({Iu\mt\Iu)) + ({Iu>\m t \u>I}) 

< ((Iu 2 \m t \Iuj2)) + ((Ioj 2 \m t \u2l}}, 

where \u>2)) is any degree-2 U(2) invariant. Since 
((Iuj2\U m ^\Iuj2)) = \{x 2 + y 2 + z 2 ) (and the equivalent 
expression for ({Iw2\U® t ' t \uj2l}}) holds for every degree- 
2 U(2) invariant, and the expression does not depend 
on t, it follows that the maximum of ({Iu>\mt\Iu))) + 
((I(jj\m t \ujl)) over all \u>)) such that ((ui\I)) = is given 
by ((Iuj2\m,2\Iu)2)) + ((Iui 2 \m2\uj2l)) . This concludes the 
proof. □ 



